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NON-BACKTRACKING LOOP SOUPS AND STATISTICAL 
MECHANICS ON SPIN NETWORKS 

FEDERICO CAMIA AND MARCIN LIS 


Abstract. We introduce and study a Markov field on the edges of a 
graph Q in dimension d > 2 whose configurations are spin networks. 
The field arises naturally as the edge-occupation field of a Poissonian 
model (a soup) of non-backtracking loops and walks characterized by 
a spatial Markov property such that, conditionally on the value of the 
edge-occupation field on a boundary set, the distributions of the loops 
and arcs on either side of the boundary are independent of each other. 
The field has a Gibbs distribution with a Hamiltonian given by a sum 
of terms which involve only edges incident on the same vertex. Its free 
energy density and other quantities can be computed exactly, and their 
critical behavior analyzed, in any dimension. 


1. Introduction 


The free energy density is an important tool and one of the main ob¬ 
jects of study in statistical mechanics, since thermodynamic functions can 
be expressed in terms of its derivatives. As a consequence, models whose 
free energy density can be computed exactly have played a crucial role in 
the development of statistical mechanics. Perhaps the main example is the 
two-dimensional Ising model, whose free energy density was famously de¬ 
rived by Onsager. However, this situation is rare and typically restricted to 
dimensions one and two, as in the case of the Ising model. 

In this paper, we introduce and study a Markov field whose free energy 
density can be computed exactly in any dimension. The field arises naturally 
as the edge-occupation field of a stochastic model of non-backtracking loops 
and walks, but it does not seem to have been studied before. Such random 
loop models first appeared in the work of Symanzik 27 on Euclidean field 
theory, and can be used to prove “triviality” (i.e., the Gaussian nature) of 
Euclidean fields in dimension five and higher 21. A lattice version appears in 
the work of Brydges, Frohlich and Spencer |6j, who develop a random walk 
representation of spin systems, and use it to prove correlation inequalities. 
The loops that appear in the work of Brydges, Frohlich and Spencer are 
allowed to backtrack. 

A prototypical example of a statistical mechanical model whose partition 
function coincides with that of a loop model is the discrete Gaussian free 
field, whose partition function can be easily expressed as the grand-canonical 
partition function of an “ideal” lattice gas of loops, i.e., a Poissonian ensem¬ 
ble of lattice loops. This is a special case of the random walk representation 
of Brydges, Frohlich and Spencer 6 . In this case, the obvious triviality of 
the Gaussian free field is reflected in the Poissonian nature of the ensemble 
of loops. In general, the ensembles of random paths and loops that appear 
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in the work of Brydges, Frohlich and Spencer are not Poissonian: the un¬ 
derlying Poisson distribution is “tilted” by means of a potential that is a 
function of the occupation field generated by the random loops and paths 
on the vertices of the lattice. This suggests that the occupation field of a 
Poissonian ensemble of loops is an interesting object to study. Such an oc¬ 
cupation field is the main object of interest in this paper, although we look 
at edge occupation and our loops have a non-backtracking condition that is 
not present in the work of Brydges, Frohlich and Spencer. 

The Poissonian ensemble of loops implicit in the work of Symanzik was 
rediscovered by Lawler and Werner 17 in connection with conformal invari¬ 


ance and the Schramm-Loewner Evolution (SLE), and named Brownian loop 
soup. A discrete version, called random walk loop soup was introduced by 
Lawler and Trujillo Ferreras jl6] , and studied extensively by Le Jan 18,19 


who discovered a connection between the vertex-occupation field of the loop 
soup and the square of the discrete Gaussian free field. For a particular 
value of the intensity of the Poisson process, the random walk loop soup is 
exactly the lattice loop model mentioned above, whose partition function 
coincides with that of the discrete Gaussian free field (up to a multiplicative 
constant). 

In this paper, we introduce a special ingredient to the the loop soup recipe: 
a non-backtracking condition on the loops. It turns out that the partition 
function of the non-backtracking loop soup can still be linked to that of the 
Gaussian free field, but the link is less direct in this case, and the proof uses 
a connection with the Ihara (edge) zeta function |12,14j[26,28j, as explained 
in Section [H 

We show that the non-backtracking loop soup possesses a spatial Markov 
property , discussed in the next sectioiQ As a consequence of this property, 
the loop soup induces a Markovian occupation field on the edges of the 
lattice or graph where it is defined, which in turn implies that the field has 
a Gibbs distribution. Quite remarkably, the Hamiltonian can be written 
explicitly and appears to be a sum of local terms, and the partition function 
and the free energy density of the occupation field can all be computed 
exactly in any dimension. These computations rely on the fact that, for 
this as for any loop soup, the partition function can be written in terms of 
a determinant. In the translation invariant case, the relevant determinant 
can be calculated, and the free energy density of the field written in closed 
form, on the d-dimensional torus for any d, as mentioned at the beginning 
of this introduction. 

Although with a different diffusion constant, the random non-backtracking 
walk satisfies a functional central limit theorem like the simple random walk 
1111. Thus, in two dimensions, in the scaling limit, the non-backtracking loop 
model should be related to the Brownian loop soup of Lawler and Werner, 
and thus to SLE. However, in view of the role that exactly-solvable models 
have played in one and two dimensions, the model we introduce in this paper 


1 Reference 


29 


posted shortly after this paper, shows that the non-backtracking con¬ 
dition is not necessary to have the spatial Markov property. The Markovian nature of 
the occupation field of t he loop soup without the non-backtracking condition was proved 
independently in 19,21 . 
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may be even more interesting in higher dimensions, since it may offer some 
insight into critical behavior in three and four dimensions, where conformal 
invariance does not play the same role as in two dimensions. 

It is also worth mentioning the intriguing connection with spin networks , 
a concept introduced by R. Penrose [24] which has applications to quantum 
gravity [251 and appears also in conformal and topological quantum field 
theory (see, e.g., [9,10])- The non-backtracking loop soup studied in this 
paper generates spin network configurations with a certain Gibbs distribu¬ 
tion, so it can indeed be regarded as a statistical mechanical model on spin 
networks. 


2. Discussion of the main results 


We introduce and study a Poissonian model of non-backtracking loops 
and arcs on graphs in dimension d > 2, and the associated edge-occupation 
field, given by the number of visits to each edge by the collection of loops 
and arcs. Adopting the language of |[lT] and [l5], we refer to this model 
as a loop soup, even though the “soup” contains in general both loops and 
arcs that start and end at boundary edges. For simplicity, in this section 
we focus only on loops; precise and more general definitions are given in the 
next section. 

Consider a connected simple graph Q = (V, E) and associate a positive 
weight x e to each edge e £ E. A (non-backtracking) loop l is a closed walk 
on the edges of Q considered up to cyclic shifts and reversal. A loop is said 
to have multiplicity me if it can be obtained as the concatenation of me 
copies of the same loop, and me is the largest such number. A loop t with 
multiplicity me is assigned weight p(£) = Y\ x e , where the product is 
taken over all edges traversed by A The loop soup £ studied in this paper is 
a Poisson point process on the space of non-backtracking loops with intensity 
measure p. The edge-occupation field Np = (Ns(e)) e£ E induced by the loop 
soup is given by the total number of visits of the loops to each edge of Q. 

Our main results concern both the loop soup itself and the induced edge- 
occupation field. 


(1) 

( 2 ) 


The loop soup has a spatial Markov property such that, conditionally 
on the value of the edge-occupation field on a boundary set, the 
distributions of the loops and arcs on either side of the boundary 
are independent of each other (Theorem 3.3). 


The edge-occupation field has a Gibbs distribution (Theorem 3.1) 
whose Hamiltonian can be written explicitly (Equation (3.1)) and 
whose partition function can be expressed as 


z = E n g n 

N vGV e£E 


N{e) 

Te 

N{e)\ 


where the sum runs over all spin network |24j configurations N and 
C v = C V (N) is a suitable function of ( N(e)) eBv . 

(3) The partition function can be expressed as a determinant and related 
to the Ihara zeta function and to the partition function of a discrete 
Gaussian free field (Section [4]). 
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(4) In the homogeneous case, x e = x, on a finite d-regular graph, Z < oo 
if x < l/(d — 1) and Z diverges for x = 1 /{d — 1) (Corollary |4.4[ ). 

In the case of translation invariant weights x e on the d-dimensional torus, 
we have the following exact results. 

(5) The partition function of the occupation held can be computed ex¬ 
plicitly in any dimension d > 2 (Corollary |5.2[ ). 

(6) The free energy density / of the occupation held on Z d can be com¬ 
puted explicitly for any dimension d > 2 as the thermodynamic limit 
of the free energy density on the d-dimensional torus (Corollary |5.3[ ). 

In the homogeneous case, x e = x Ve E E, we have the following exact results 
on Z rf . 

(7) The d -dimensional free energy density f(x) is hnite for x < 1/(2 d — 
1), and has a singularity as x /* 1/(2 d — 1) which depends on the 
dimension d (Corollary |5.4[ ). 

(8) The truncated two-point function decays exponentially for x < 1/(2 d— 
1) (Corollary |7.3[ ). 

Points (7) and (8) show that the edge-occupation held undergoes a sharp 
phase transition and that the critical point can be explicitly computed and 
the critical behavior analyzed for periodic graphs and lattices. This is done 
by studying the free energy density in the thermodynamic limit and by 
deriving expressions for the truncated two-point function. 

In addition to the results mentioned above, Sections [6] and [7] contain more 
results on the distribution and on the two-point function of the occupation 
held, respectively. 

In the last section of the paper, we use the loop soup to dehne a spin 
model on the vertices of the dual graph. The spin model can be shown to be 
reflection positive using to the Markov nature of the edge-occupation held. 
In two dimensions, we conjecture that the scaling limit of the spin held is 
one of the conformal helds introduced in |T] (see also j8j). We note that the 
same spin model generated by an ordinary loop soup was introduced by Le 
Jan [l9| and is also reflection-positive. 

3. The non-backtracking loop soup 

Let Q = (17, E) be a connected simple graph, and let E be the set of its 
directed edges. For a directed edge e = {tg, hg ) E E, —e = ( hg, tg ) E E is its 
reversal, and e = {tg, hg} E E its undirected version. We assume that the 
graph is equipped with a positive edge weight x e for each e E E. By dQ we 
denote the boundary of Q , i.e. the (possibly empty) set of edges incident on 
a vertex of degree one. 

A (non-backtracking) walk u of length |w| = n > 1 is a sequence of 
directed edges u = (coi,..., w n +i) such that Lj i+1 = h^ i and uJj+i / — Wj 
for 1 < i < n. Note that the length of a walk is the number of steps the 
walk makes between the edges rather than the number of edges itself. By 
cJ -1 = (—cJ|£|+i,..., — cui) we denote its reversal, and for two walks uj, Co' 
such that W|£| +1 = uj{, we define the concatenation 

IX © Co' = (wi, . . • , W2) • ■ ■ ) dj\ cJ'l+l)- 
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Rooted, loops are walks starting and ending at the same directed edge. 
The multiplicity of a rooted loop Co, denoted by rn$, is the largest number 
m such that Co is the m-fold concatenation of some rooted loop c o' with itself. 
An unoriented walk, is a walk without a specified direction of traversal, i.e. 
a two-element equivalence class under the relation Co ~ cW 1 . An arc is 
an unoriented walk starting and ending on dG. Arcs will be denoted by a. 
Unrooted loops are equivalence classes of loops under the cyclic shift relation 
Co ~ (ui, iWi+i, ..., W|£|+i, Co\, C 02 , ■ ■ ■ cJi- 1 ). Unrooted unoriented loops will be 
simply referred to as loops and will be denoted by £. 

With a slight abuse of notation, if a function / defined on walks is in¬ 
variant under reversal, then f(a) is the evaluation of / at any of the two 
representatives of the arc a. Similarly, if a function g defined on rooted 
loops is invariant under reversal and cyclic shift, g(£) is the evaluation of g 
at any representative of l. The weight of a walk Co is 


|3| 

*(w)= n v x ^ x “i+i- 
2=1 


The loop and arc measures are given by 

x(£) 


g(£) = 


and 


me 


Udi.cn) 


x(a). 


By L we will denote a realization of a Poisson point process with intensity 
measure fi, and by A a realization of a Poisson point process with intensity 
measure (jlq. We will write S = £U A, where C and A are independent. The 
partition function of C is 


z c =j2 w i L )=j2 n 

l l eeL 


g {£)* 1 

Wv 1 


where the sum is taken over all multi-sets L of loops, called loop configura¬ 
tions, and where #£ is the number of occurrences of £ in L. Similarly, the 
partition function of A is 


z A =Y, w i A )= y. n 




Ud{a)* a 

(#«)■ 


where the sum is taken over all multi-sets of arcs A, called arc configurations. 
The partition function of S is given by Z$ = ZnZ A . We will only consider 
the cases where Z^Z A < oo. A multi-set of loops and arcs will be called 
a soup configuration or simply a configuration. In particular, loop and arc 
configurations are soup configurations. 

For a soup configuration S, we define Ns = ( Ns{e)) e& E to be the network 
(or edge-occupation field ) induced by S, i.e. the total number of visits of the 
walks from S to each edge of Q. Ng is a spin network in the sense of Penrose 
24 , that is, if the edges incident on vertex v are e\,... ,ek, then Ya=i ^s{ci) 
must be even and that it cannot be smaller than 2maxj = x i ... ) fe Ns(ei). 

One of the main results of this paper is that the distribution of the ran¬ 
dom network N$ with prescribed boundary conditions £ = (A r (e)) eS ag (in 
particular, the distribution of Nc for zero boundary condition) is given by 
a Gibbs distribution with a local Hamiltonian. 
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More precisely, suppose that e\,... ,ek are all the edges incident on vertex 
v and imagine replacing edge e* by Ng{ei ) distinct, colored, edges. Assume 
that each colored edge incident on v has a unique color and let C v be the 
number of different ways in which those edges can be connected in such a 
way that each colored edge corresponding to e* is connected to a colored 
edge corresponding to ej for some j / i. (It is clear that one can always 
connect all colored edges in this way since, after all, Ng(ei) is the number 
of visits to edge e* of the loops and arcs from the soup.) If Q = (V. E) is 
finite, we can define the Hamiltonian H = H{Ng) = YIvgv w h ere ) for 
each vertex v with edges e\, ..., e*, incident on it, 


(3.1) 


H v = - log 


c v 


=n 

/ rii=i IVs(ej)! *=' 


,N s {ei)/2 


Theorem 3.1 (Gibbs distribution). If the graph Q is finite and S is a soup 
of loops and arcs in Q such that Zg < oo, then the distribution induced 
on the edge-occupation field by S is the same as a Gibbs distribution with 
Hamiltonian H. That is, the edge-occupation configuration Ng has probabil¬ 
ity 

1 -H{N S ) 

Z 

where 

N(e) 

(3.2) Z = £ e~ Hm = E II C » II my = 

N N vev e&E v '' 

with the sums running over all network configurations N. 


We note that the partition function (3.2) is reminiscent of the random 
current representation used by Aizenman m (see also |3] for a more recent 
application). 

If Q is a trivalent graph, the combinatorics is simple and one can easily 
compute C v . (R. Penrose uses trivalent graphs to define spin networks pre¬ 
cisely for this reason, although the concept is more general-see [ 241.) Let 
ei, 62 , c-i be the edges incident on some vertex v of Q, with occupation val¬ 
ues N(e i), N(e 2 ) and JV(e 3 ), and define N 12 = N ^ ei '> +N ^ 2 '>~ N ^ e3 ) ; ]\j 13 = 
Af(ei)+Af(e 3 )-IV(e 2 ) ^ _ N(e 2 )+N(e 3 )-N(e 1 ) ^ are ^g un iq u g solutions 

of the system of linear equations ]Vi 2 + A 13 = iV(ei), N± 2 + -/V 23 = N(e 2 ), 
N 13 + N 23 = N(e 3 ).) It is easy to see that C v = and 

consequently 


H v = - log 


( V N ( e lV- N (e2) ] -N( e 3)l N(ei )/2 N(e 2 )/2 N(e a )/2\ 

1 lL e2 ^es I 


[ N 12 !N 13 <N 23 ! 


Theorem 3.1 is a consequence of the conditional factorization property 
of the law of S, which implies that the soup possesses a spatial Markov 
property. In other words, the edge-occupation field is a Markov field. A 
natural way to express this property is to cut loops into arcs, and we will 
now make precise what we mean by that. Suppose that we fix a set of edges 
H C E. By Qh w e denote the modified graph Q where the edges from H 
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Figure 1. Graphs Q with H = {ei, 62 , 63 }, Qh, and Q' H for 
S with N s {e 1 ) = 2, N s (e 2 ) = 0, N s (e 1 ) = 2 


are cut in half, i.e. for each e = {u, v} E H, we remove e from the edge-set, 
and add two new edges, called half-edges , {u,u'}, {?/,•?;} with v! 7 ^ v' , and 
i/, r/ ^ V. The weight of each half-edge is equal to the weight of the removed 
edge it replaces. Note that the half-edges belong to OQh . Each loop or arc 
in Q visiting H , when cut along H , gives rise to a multi-set of arcs in Qh, 
corresponding to its maximal subwalks which visit only edges from E\ H, 
with the possible exception at the endpoints. (These arcs are excursions the 
walks make outside of H.) If a walk does not visit H , then it is not affected 
when the edges of H are cut in two. Given a configuration S in Q, we define 
Sg H to be the configuration in Qh resulting from cutting the walks from S 
(taken with multiplicities) along H. 

Lemma 3.2 (Conditional factorization). For a set of edges H , and a con¬ 
figuration S, 

Y w(T) = w(Sg H ) JJ N s (e)\ 

T ■ T e H =Sg H eeH 

Proof. Note that all T with Tg H = Sg H differ only in the way the arcs of Tg H 
are connected to each other along the edges in H. It is enough to consider 
the case when all walks in S visit H and hence Sg H consists only of arcs. 
This is because the contribution of the loops from T n Tg H to w (T ) is the 
same for all T with Tg H = Sg H and is equal to the total loop contribution 
to w(Sg H ). 

We will prove the result by showing that both the l.h.s. and the r.h.s. of 
the equality are proportional to the number of ways one can assign colors 
to the visits of walks of a configuration to the edges of H. In order to count 
that number, let Q’ H be the extended graph Qh- where for each e E H, we 
connect the degree-one vertices of the two corresponding half-edges in Qh 
with Ns(e) distinguishable parallel edges (see Figure [T]). We can think of 
the new edges as being colored by Ns(e) different colors. 

Consider configurations in Q' H which visit each of the colored edges exactly 
once. Each such configuration maps to a configuration in Q by forgetting 
the colors, i.e. identifying the colored edges for each e E H. Note that the 
resulting configuration T in Q satisfies Nt\h = Ns\h- This mapping is 
many-to-one, and the cardinality of the preimage of each configuration T 
in G is proportional to w(T). We will now justify this statement. Fix a 
configuration T in Q such that Tg H = Sg H . Consider a set T where each 
repeated loop or arc in T is made distinguishable by adding some additional 
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markers, and moreover where for each loop we choose a root and an orien¬ 
tation. In this case, since all steps of the walks in T are distinguishable, we 
have Ns(e)l different ways of assigning colors to the visits of the walks in T 
to e € H. What is left to do is to account for the multiplicities coming from 
identifying the marked loops corresponding to the same unmarked loop and 
the rooted directed loops corresponding to the same unrooted undirected 
loop. This gives 


,, „ rwA-sM! = iwvgW 1 „. m 

IUt(#«)!IW(#«)! mf IW IL CT aK°0 # “ 

different colorings of the visits to H of the walks in the configuration T. 
Note that here we have used the fact that each non-backtracking walk has 
two distinct oriented versions. 

Thus, the number of ways one can assign colors to the visits of walks of 
a configuration to the edges of H can be written as 


(3.4) 


E 

T: Tg H =Sg H 


_ w{T) _ 

n« e r *(*)** Ilaer *(“)*' 


n N s( e v- 

eeH 


Let us now derive a different expression for this number, this time by 
counting the number of ways one can connect the arcs from Sg H into loops 
in Q' h in such a way that each colored edge is used only once. To this end, 
take an edge e = {u, v} from H and consider the corresponding two directed 
half-edges e\ = (v',v), e ?2 = (u',u), and Ns(e ) colored edges in Q' H . For the 
purpose of this proof we will consider directed arcs. Again, the property of 
non-backtracking walks that we use is that each undirected walk has two 
distinct directed versions. Let Sg H be the multi-set of all directed versions 
of the arcs from Sg H , and let A(ei) be the multi-set of directed arcs from 
Sg H starting at e*, i = 1, 2. One has 


^2 #® = N s( e ) for 2 = 1,2, 

a£A(ei) 


where #a is the multiplicity of the directed arc a in A(e*) (which is equal 
to the multiplicity of its undirected version a in Sg H ). We now distribute 
the Ns(e ) colors between the directed arcs in A(ej), * = 1,2. Since the arcs 
have multiplicities, there are exactly 

N s (e)\ x jVg(e)! 
na ie T(ei)(#“l) ! n 5 2e A(e2)(^“ 2 )' 

such assignments. If we take the product over H and use the fact that each 
arc has two directed versions, we arrive at 

n, 6H iVs(4)! 2 

where is the multiplicity of a in Sg H . This is the number of all possi¬ 
ble assignments of colors to the directed arcs. We now want to forget the 
orientation of the arcs, so, for each arc a, we need to pair up the opposite 
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directed, colored arcs a. and a 1 . Since we can pair any colored arc a with 
any colored arc a -1 , we have {ffa)\ different pairings. Hence, 

11 ez H Ns(e)\ 2 _ w(Sg H ) 


(3.5) 


n c 


La£Sg H ~ v “/ e&H 

is the number of all possible ways of connecting the arcs from Sh in such a 
way that each colored edge is used once, which gives us another expression 
for (3.4). 

We remind the reader that we need only consider the case when all walks 
in S visit H and hence Sg H consists only of arcs, and note that, for every T 
such that Tg H = Sg H , we have the identity 

Yl x(a)* a . 

a£Sg H 


x(a)# c 




n^) #£ n 

feT aeT 


ar = 


Using this identity and comparing (3.4) with (3.5) concludes the proof of 
the lemma. □ 

Given a subgraph Q\ of Qh which is a union of a number of connected 
components of Qh , we define Sg 1 to be Sg H restricted to walks in Q\ . If £ : 
dQ — > IN>o are boundary conditions, then we write Pg £ for the probability 
measure governing S defined on G and conditioned to satisfy N$\gg = £. 

Theorem 3.3 (Spatial Markov property). Let H be a subset of edges of Q, 
and let Q\ be one of the connected components of Qh- Let S be a configura¬ 
tion in Q. Then, for all boundary conditions £ = (N(e)) e& gg on dQ, 

(‘-’di ■~4?i | Sg H \gi — *S , g ff \5 1 ) 

= Pg,£ (Sg t = S'eJ A(s|aei = ^sladi) 

= P g 1 ,s 1 (S = Sg 1 ), 

where £i are boundary conditions on dQ\ given by 

£ on dQifldQ, 


= 


NsIdSi 071 dQi \ dQ. 


Proof. From Lemma 3.2 and the factorization property of the Poisson point 
process weights, it follows that 

= S Bl \ Sg H \ Gl = Sg H \ Gl ) 

= P G ,{ {Sg H = Sg H )/ P G ,f = S Gh\Gi) 

= w ( T )/ E u, ( T ) 

T'-Tg H —Sg H T.Tg H \g 1 —Sg H \g 1 

= w{S Gl U Sg H \ gi )/ W ( T 1 U 

Ti:N Tl \dg 1 =Ns\og 1 

= w{S Ql )/ Y, w{T 1 ) = Pg Ul {S = Sg 1 ), 

Ti:N Tl |asi = N s\ag 1 

where T\ denotes a configuration in Q\. □ 


We conclude this section with a proof of Theorem 3.1 
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Proof of Theorem \3.1\ Let N = {N(e)} e £E be a collection of numbers in 
!N>o indexed by the edge set E of Q. We write S ~ IV if the soup con¬ 
figuration S induces the occupation field N (i.e., if Ng = N) and define 
S(N) = {5 : S ~ IV} = (S' : Ng = IV}. Let Q' be the extended graph Q 
where each e G E is replaced with Ng(e ) distinguishable parallel edges. We 
can think of the new edges as being colored by Ng(e ) different colors. Con¬ 
sider colored configurations S* in Q' which visit each of the colored edges 
exactly once. We write S*(N ) for the set of such colored configurations and 
|5*(./V)| for the cardinality of S*(N). 

With this notation, the partition function of the soup can be written as 


Z S 



£n 

A aeA 


it 


n 

aeS 


x(a ) 

(#«) ! 


ho(a)* a \ 

(#«)! ) 


e n x e {e) e n 

N e&E S~Nt£S 


1 


n 

aes 


1 

(#«) ! 


Ei^win 

N e&E 


N(e) 


N(e)\ 


En^n 

N v£V e£E 


X 


N(e) 

e 


N(e)\ 


= En^-=E«- ww =^ 

N vev N 

where we have used the obvious identity |<S*(iV)| = Y\ v(Z yC v and the non¬ 
trivial identity 


£n 

s~n ees 


1 


n 

aes 


1 

(#a) ! 


S*(N)\ n 

e£E 


1 

IV(e)! ’ 


The latter identity can be derived using arguments similar to those in the 
discussion preceding (3.3), as explained in the lemma below. 

x(£) 


To conclude the proof, we write 


£ p s(s) = ^-eii 


S~N 


Z S S~N£gS a&S 

1 


n x(a) 
(# a )- 




e&E 


N(e) 
N(e)\ 

i N{e) 

z n n ]%)i 

v&V e&E v ' 

z 


In 


vev 


-h v 


Lemma 3.4. With the notation introduced in the proof of Theorem 3.1, the 
following identity holds: 

1 TT 1 


\s*m=u N ^y- E n 


e£E 


S~N fgs ces (# a )' 


n 
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Proof. We will prove the result by counting the number of ways one can 
assign colors to the visits of walks of a configuration to the edges of Q. To 
this end, let again Q' be the extended graph Q where each e £ E is replaced 
with Ns(e) distinguishable parallel edges. We can think of the new edges as 
being colored by N$(e) different colors. 

Consider configurations in Q' which visit each of the colored edges exactly 
once. Each such colored configuration S* maps to a configuration S in Q 
by forgetting the colors, i.e. identifying the colored edges for each e £ E. 
Note that the resulting configuration S in Q satisfies Ng = N, and that 
the mapping is many-to-one. We want to determine the cardinality of the 
preimage of each configuration S in G. 

Fix a configuration S in Q such that Ns = N. Consider a set S where each 
repeated loop or arc in S is made distinguishable by adding some additional 
markers, and moreover where for each loop we choose a root and an orien¬ 
tation. In this case, since all steps of the walks in S are distinguishable, we 
have N$(e)\ different ways of assigning colors to the visits of the walks in S 
to e £ E. What is left to do is to account for the multiplicities coming from 
identifying the marked loops corresponding to the same unmarked loop and 
the rooted directed loops corresponding to the same unrooted undirected 
loop. This gives 

(3.6) n. £E ^(e)! 

n„ £S (#a)! n K swM l 

for the cardinality of the preimage of S. Here we again used that each 
non-backtracking loop has two distinct oriented versions. (Note that the 
counting doesn’t work for ordinary loops without the non-backtracking con¬ 
dition because in that case one can have loops that have a single oriented 
version instead of two-think of a loop consisting of a single edge traversed 
twice.) Summing over all configurations S £ S(N ) concludes the proof. □ 

4. Determinantal formulas for the partition function 

In this section we express the partition function of our model in terms of 
determinants of two different matrices, which, for certain values of the edge 
weights, involve transition matrices of some Markov processes. The process 
involved in the first determinantal formula is an asymmetric random walk 
on the directed edges of Q, and the one involved in the second formula 
is a random walk on the vertices of Q. As a consequence of the second 
determinantal formula, we derive a relation between the partition function of 
our model and that of the discrete Gaussian free field. Along the way we also 
uncover a connection between the partition function of our model and the 
Ihara zeta function. In the next section we will use the first determinantal 
formulas to do exact computations. 

We assume that Q is finite and connected and has no vertex of degree 1. 
For e,g £ E, let 


(4.1) 


^-e,g ~ 


x e if hg = t(j and tg / hg, 
0 otherwise, 
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and let p( A) be the spectral radius of A. The next observation, a well known 
result, is crucial and allows for an exact solution of our model. 


Lemma 4.1. The partition function Zc is finite if and only if p{ A) < 1, in 
which case 

Zc = det - 2 (Id — A), 


where Id is the identity matrix indexed by E. 

Proof. Let A j, i = 1,2, be the eigenvalues of A, and let W5- be 

the set of all oriented walks of length n starting and ending at the oriented 
edge e. Using the definition of the loop measure, we have that for each n, 


= 

t\e\=n 


£ 

l:\i\=n 


x(i) 

me 


ST' ST' x (£) 

Z 9 79 

geSfeW2,g 


Tr A n 
2 n 


\E\ 

£ 

2=1 


2 n 


It follows that the first expression is summable over n if and only if maxj | Aj 
p(A) < 1. From the definition of the partition function, we have that 

( OO 

n= 1 i-.\t\=n 

and hence the first part of the lemma follows. To finish the proof, we notice 
that 

oo \E\ \E\ \E\ 

—2 In Zc = - EE£T = E ln(1-Ai) = ln JJ(1— Xi) = lndet(Id—A). □ 

n= 1 i= 1 i= 1 i=l 


In the next section we will use Lemma [4.1| to perform exact computations. 
First, however, we describe an interesting identity providing a different de- 
terminantal representation for the partition function. Such a determinantal 
representation was introduced in connection with the Ihara (edge) zeta func¬ 
tion |12l|14 26 


which is defined as the infinite product 


c( x ) = n a-*<*)) 


-1 




where x = (x e ) e es denotes the vector of edge weights and V is the collection 
of all oriented, unrooted, non-backtracking loops with multiplicity 1. Note 
that in the general definition of the Ihara zeta function, the weights x e can 
be complex and can be different for the two opposite orientations of an 


undirected edge. By grouping together loops in (4.2) that are multiples of 
the sa me lo op Co € V, we easily obtain that Zc = C^( x )) an d hence by 
Lemma 4.1 £ _1 (x) = det (Id — A). (This is the content of Theorem 3 of 
and Theorem 3.3 of fl3 


26 


To state the alternative determinantal formula, we need to first define 
two matrices indexed by the vertices of Q. If Hx^ < 1, let D be a diagonal 
matrix with entries 
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and let 


Av,u — 


i-®! 


if e = {u, v} G E, 
otherwise 


be a weighted adjacency matrix. 

Lemma 4.2. For HxHoo small enough, 

Z c = det _1/2 (Id + D - A) JJ (1 

e&E 

where Id is the identity matrix indexed by V. 


- x, 


2 \— 1/2 


Proof. This result follows from Theorem 2 of 128 and Lemma 4.1 


□ 


We will use Lemma 4.2 to relate the partition function of our model to 
that of a discrete Gaussian free field. To this end, attach a weight c e > 0 
to each edge e G E and a killing rate k v > 0 to each vertex v G V, and let 
A„ = YleBv c e + k v . The weights and killing rates induce a sub-Markovian 
transition matrix P on the vertices of Q with transition probabilities p VtU = 
if e = {f, u} G E, and 0 otherwise. The transition matrix is A-symmetric, 
meaning that A v p VjU = A u Pu,v One can introduce a “cemetery” state A and 
extend P to a Markovian transition matrix P on kUA by setting p v a = 
and pa,A = 1- We assume that the Green’s function Gp(v, u ) of the Markov 
chain associated to P (i.e., the expected number of visits to u of the chain 
started at v) is finite. This condition is equivalent to saying that p(P) < 1 
or that the Markov chain is always absorbed in A. 

The discrete Gaussian free field on Q associated with the transition ma¬ 
trix P is a collection {<f v ) v ^v °f mean-zero Gaussian random variables with 
covariance E ((j) v 4> u ) = Gp(v,u). It has a Gibbs distribution e~ H p /Zp FF 
with Hamiltonian 


ztGFF 

rLp 


(p) = l E 

{v,u}eE 


1 


^{v,u}{^V U ) cy ^ ^ 1 


vev 


from which it follows that its partition function, Zp F F , is given by a Gauss¬ 
ian integral and can therefore be represented in terms of a determinant, as 
follows: 


(4.3) 


Z' 


GFF 


n 

vev 


( 27T\1/2 1 

' ' det ~2 (Id — P) 


vaJ 


where Id is the identity matrix indexed by V. 

Theorem 4.3. Let x be a vector of edge weights such that HxH^ < 1 and 


(4.4) 

For v & V, let 

Then, the matrix 


E x e 
1 + X t 

e^v 


<1 \/v G V. 


P* = 


A* = 1 


_1_ Xe 

AJ 1—2 

0 


+ £r 


e3v 


~ Xt 


if e = {v, u} G E, 
otherwise 
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is a X* -symmetric, sub-Markovian transition matrix on V. Moreover, if 
Z c < oo, 


(4.5) 


Z c = (2vr)-l v l/ 2 Ha- x l)~ 1/2 Z P* F i 


e£E 


where P* is the Markovian extension of P*. 

Proof. The matrix P* is obviously A*-symmetric and an easy computation 
shows that the inequalities 


e3v 


+ X e 


< 1 


and 


E l x, 

X* 1 

e3v v 


— X 


2 — 


< 1 


are equivalent. Hence, condition (4.4) guarantees that P* is sub-Markovian. 


Using Lemma |4.2[ for ||x||oo sufficiently small, we obtain 

{Zc)- 2 = n (! - x e) det ( Id + D - A) 

„2 


e&E 


n<i - 4) n i 1 +S rv^) det < Id - p ”> ■ 


e£E 


v£V 


e3v 


By Lemma 4.1 the l.h.s. and the r.h.s. are both polynomials, so they are 
equal for all x. Hence, 

Zc = IL 1 - *2F 1/2 n(‘ + E ^("^detAdd - P*) 

eOE veV e3v ^ 

for all x such that 0 < Zc < oo. Using the definition of X* and the deter- 
minantal formula 4.3 concludes the proof. 

□ 

Corollary 4.4. Let x be a vector of edge weights such that HxHoo < 1 and 

x e 


S r 


e3v 


+ X e 


IfBv^V such that 


S T 

e3v 


< 1 Vv E V. 


< 1 , 


+ x t 


then Zc < oo and Equation ( |4.5[ ) holds. If 

x 


Sr 

e3v 


+ X e 


= 1 Vu £ V, 


then Zc diverges. 

In particular, in the homogeneous case, x e = x, on a d-regular graph, 
Zc < oo if x < 1 /{d— 1) and Zc diverges if x = 1 /{d— 1). 

Proof. It is easy to verify that, if YheBv T+f~ — ^ f° r ad v £ V and moreover 
Yle 3 v J+P~ < 1 for at least one v € V, then Zc < oo and p(P*) < 1. If, 
in addition, || x lloo < 1, then Theorem 4.3 and the determinantal formula 
(4.3) hold. It is also easy to check that, if YleBv UvU = ^ ^ or ad v e the 
transition matrix P* is Markovian and Zc diverges. □ 
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Remark 4.5. A direct way of proving the homogeneous case result is by 
counting all possible non-backtracking walks on the graph, or by noticing 
that p( A) < || A|j i = x(d— 1) and using Lemma 4.1 Similarly, for any graph 


with vertices of degree d v > 2 , one can show that a sufficient condition to 
ensure that Zg < oo is: x e < l/(d v — 1) for all v E V and all e 3 v. 

We conclude this section with an interesting observation. As mentioned in 
the introduction, the partition function of the Gaussian free field in the form 


(4.3) can be easily expressed in terms of the partition function of a “regular” 
random walk loop soup, i.e., without the non-backtracking condition. (The 
proof of this fact is completely analogous to the proof of Lemma 4.1 The 
interested reader is referred to the proof of Lemma 1.2 of 1 6 ].) This provides 
a link between the partition function of the non-backtracking loop soup 
with transition probabilities given by A and the random walk loop soup 
with transition probabilities given by P*. It is possible that the connection 
between those two loop soups and their occupation fields is realized at a 
deeper level than that of the partition functions. Such a deeper connection, 
if it exists, could be revealed by an analysis of the determinantal formula 
for the Ihara zeta function given in Theorem 2 of 1281. 


5. Exact computations 


In this section we will compute explicitly the free energy density of trans¬ 
lation invariant models and the one-point function of homogeneous models 
on the torus /{nZj) d for any n > 1 , and d > 1 , and hence after taking 
the thermodynamic limit, on all hypercubic lattices d > 1. We will 
also prove the exponential decay of the two-point function in the subcritical 
regime. The reason why an exact solution is available is the fact that the 
partition function of the model is given by the square root of the determinant 
of a matrix (a situation similar to the one of the Ising and dimer model, and 
also the discrete Gaussian free field). It follows that the relevant quantities 
are expressed in terms of the eigenvalues and can be computed for periodic 
graphs using the Fourier transform. Unlike in the Ising and dimer model 
case, the determinantal formulas are valid in all dimensions, similarly to the 
case of the discrete Gaussian free field. 


5.1. Partition function and free energy density. Let g* = Z d /(nZ) d 
be a d-dimensional torus of size n d . For a vertex k = (k\, & 2 , • • •, k n ) E G d , 
and a unit direction vector v = (iq, V 2 , ■ ■ ■, v n ) such that Vj = ±1 for some 
j, and Vi = 0 for i ^ j, we will write (k, v) for the directed edge e with 
tg = k and hg = k + v. We consider a translation invariant weight vector 
x = (aq,..., xg) that assigns weight Xj to undirected versions of all edges of 
the form (k, v), where v is a unit vector in the j-th direction. In this case, 


the transition matrix (4.1) is given by 


(5.1) 


Kit 


k.v),(l,u) 


0 


if k + v = 1 and v ^ — u, 
otherwise. 


Consider its Fourier transform 

[\J(p,v),(q,u) = ^2 ' 

k,ie^ 


3si(p-k—q-l)r A d 


IKk 


k,v),(l,u) • 
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It is easily seen that 

[An](p,v),(q , u ) 


_ ) 3 


Xje n 


p v if p = q and v 7 ^ —u, 
otherwise. 


Note that A^(p) is block-diagonal with blocks indexed by the vertices p £ 
whose rows and columns correspond to the directed edges e satisfying tg= p. 
Let A£(p) denote the 2d x 2d block corresponding to p € Since A.% is 
similar to A^. one has that 

(5.2) det(Id — A^) = det(Id — A^) = det(Id — A^(p)), 

where Id is the 2d-dimensional identity. One can explicitly compute these 
determinants, as shown below. 


Lemma 5.1. We have that 

1 + 2 ^^(^- cos (f^)) n 


det(Id — A^(p)) = 


i =1 


1 — X- 


= P. n . 


na-f)- 

i= 1 

For d = 1, 


Id - A* (p) = 


Proof. The proof is by induction on d. Let Zj = e 

1 — x\Z\ 0 
0 1 — x\Z\ 

and the statement is true. Assume that it holds true for d < k, and consider 
the matrix A^ + 1 (p) for p € G% +1 - Let p' = (pi,P 2 , ■ ■ ■ ,Pk) £ Gn be the 
restriction of p to the first k coordinates. For a number a, we will write 
a for a row- or column-vector with entries all equal to a. Let M^(p) be 
the matrix Id — A^(p) where each row corresponding to a directed edge 
(p, v) is divided by —.(•+■ « p v . where p • v is either pj or —pj. Note that 
det(Id — A^(p)) = det M% (p) nii xf. Using linearity of the determinant, 


we have 


detAf* + 1 (p) = 


•J^ _ z k -\-1 

x k -\-1 

0 

1 


0 

y _ z k -\-1 

x k+1 

1 


1 

1 

M*(p') 


y _ z k-\-l 

x k -\-1 

0 

1 


-1 


= 1 - 


= 1 - 


+ 1 - 


= 2 - 


Zk+1 

Xk+1 

Zk+1 

Xk+1 

Zk+1 
Xk+1 
2 

Xk+1 


z k -\-1 
_ x k +1 

y _ z k-\-l 

x k -\-1 

1 

y z k-\-1 

x k +1 

l 

l 1 

I M k {p>) 


0 

1 

M^p’] 

1 

M%( p') 


+ 1 - 


Zk+l \ 


%k~\~ 1 




+ 


1 _ Zk + i \ 


_ z k -\-1 
x k -\-1 

1 


0 I 

1 M*( p') 

0 

m*( p ') 


1 

M n(p') 

2iri 


+ 1 - 


cos 


/27TZ 

V n 


Pk+l 


Zk+1 \ 
Xk+1 ' 

1 

*n (P' 


AT*'-'' 


-1 

1 

+ 


p'; 


J k+1 


- 1) detM*(p'), 
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Let Mn(p') be the matrix M k (jp'), where from each entry we subtract 1. It 
is a block diagonal matrix with blocks of size 2 : 


Xi 

_1 _£i I ’ 


for 1 < i < k, whose rows and columns correspond to the pairs of directed 
edges (p',±v). Hence, 


1 


1 1 


i i 

I m*(p') 


0 M*( p') 


detM^(p / ) = n(L 2 - 1 

2=1 L 


Therefore, by the induction assumption, 
det(Id - A£ + 1 (p)) 

k 

r / ’/'ir \ ~\ . 

= 2 xfc + i 


2vr 

•^fc + 1 ^OS l Pfc + 1 
n 


IK 1 - ) + (1 - 4+l) det ( Id - ^n(P)) 


2=1 


/c+1 


! + 2 E 


- COS (yfPi)) 


2=1 


1 - xi 


k+1 


mw)- 


□ 


2=1 


From all these considerations we obtain an exact formula for the partition 
function of the model on the torus. 

Corollary 5.2. The partition function of the model on with translation 
invariant weights x = (xi ,..., Xd) satisfying 

(5-3) f 


2=1 


1 

1 + X{ 2 


is 


z c=n 

pee;! L 


d ( /2ir V)1 2 fil 

l + 008 (tT**)) 


2=1 


l — xf 


IK 1 - %i) 


r, n d 
Z ' 2 


i=l 


Proof. We use (5.1) and Lemma 


blocks are positive whenever (5.3) holds true. 


5.1[ and note that the determinants of all 

□ 


The free energy density of the model is defined as minus the logarithm of 
the partition function divided by the “volume” (the number of edges): 

log Z c 


/(*) = - 


\E\ 


As an easy consequence of the corollary above, we obtain that the limiting 
free energy density as Qf, approaches is given by an explicit formula. 

Corollary 5.3. The free energy density of the model on with translation 
invariant weights x satisfying (5.3) in the thermodynamic limit /* ’ZA is 
given by 

/(x) -2 d = Y^ log(l - xj ) + [ log 


2=1 


1 

W) 


'[ 0 , 271 -]° 


xAxi — COSO;) 

1 + 2 E "■ 


2=1 


l — xi 


da. 
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Note that the logarithm in the integral diverges as a -A 0 on the critical 
surface 

d 

Sr 


Xi 


i =1 


+ Xi 


From now on, we will simplify the setting by considering only homogenous 
models with a single parameter x such that x = (x ,..., x). In this case, the 
critical point is x c = 1/(2 d— 1). We now analyze the behavior of the singular 
part of the free energy density as x /"• x c . We will write A(x) ~ B(x) 
if c\B(x ) + C 2 < A(x) < C\B(x) + C* 2 , for some constants ci,C 2 ,Ci,C 2 
(depending on d) as x x c . 

Corollary 5.4. Let 5 = d/2 if d is even and 6 = (d+l)/2 if d is odd. If n < 
<5, then d n f(x)/dx n stays finite as x /*■ x c . If d is even, then d s f(x)/dx s ~ 
log(x c — x), and if d is odd, then d s f(x)/dx 5 ~ (x c — x ) _1//2 as x /* x c . 

Proof. Letting 2 xp(x) = (2d — l)x 2 — 2dx + l (so that p(x c ) = 0 and 1 
const / 0 as a; / x c ), up to constants, one can write the singular part of / 
as follows: 


'[0,27r] d 


log 


= 2 


'to, 


p(x) + Ed — cos af) da 
i =1 

d 

log p(x) + Y^(l — cos ai) 

■\ A * ^ 


i= 1 
d 


da 


log 


i— 1 


[0,7r/2] d 

r _ 

log \p(x) + ^yA 


p(x) + Ed — cos ai) da 
dy% 


= 2 


f 

'[ 0 ,l] d 

- f 


i =1 
d 


= l Vhii 2 - Vi) 


log 


p(x) + 'Yh yi n y i 1/2dyi 


i= 1 i= 1 

d d 


log p(x) + J 2 z i n dz, 


1=1 


J[ 0,1] d 
rVd 

^ ' log \p(x) + r 2 ]r d ~ 1 dr 


i= 1 


'0 


= 2 d d d / 2 ~ l log (p(x) + d\ - f 

d Jo 


r\fd r d +1 


d Jo p(x) + r 2 


dr. 


We see that the integral is convergent at x = x c for all d > 2. However, 
taking n derivatives of / generates a term containing the integral 


r'Jd 


r d 1 dr 


Jo 0 p(x) + r 2 ) n ' 

At x = x c , this integral is convergent if 2 n < d and divergent if 2 n > d. 
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Writing p = p(x), for x sufficiently close to x c , we have that 


and 


/•Vd r d-i dr 

lo (P + r2 ) n 


f 

■Wp 


fd+P (s 2 -p)- 2 ~ l 


a 2n-1 


-ds < 


j-Vd+p 

/ ' 
Wp 


„d—l—2n 


ds 


r-y/d r d~ 1 dr r d ~ 1 dr r d ~ 1 dr 


I o (P + r 2 ) n J o (p + r 2 ) n J^p (p + r 2 ) r 


> 


1 


rVP 


r d_1 dr + — 


ry/d 


(2 P) n Jo 

1 


„d— 1— 2n 


dr 


>Vp 


= -v —P 

dT 


d/2—n 


+ 


("sfd 


'Vp 


„d—l—2n 


dr. 


The last statement of the lemma follows taking n = 5. 


□ 


5.2. The one-point function. In this section we compute the one-point 
function of the homogenous model on Q d and Z d . We begin with a lemma 
which expresses it in terms of the Green’s function. The result is proved by 
expressing the desired quantity in terms of a similar object in the soup of 
oriented loops , and then repeating a classic proof of an analogous statement 
for general loop soups [l9|. Let 

Ggg = [(Id - A) -1 ]^ 

be the Green’s function for the non-backtracking random walk. If X is a 
random variable, we will write (X) for its expectation. 


Lemma 5.5. For any edge e, 

(N c (e)) = Gg t g — 1, 

where e is any of the two oriented versions of e. As a consequence, 


(w\^ Ndt 


e£E 


1 


e = 


Tr(Id - A )" 1 - 1 


2\E\ 


Proof. Let C be the soup of unrooted but oriented loops with intensity 
i.e. a Poisson point process with intensity measure Jfl, where fl(£) = p(£), 
and where l is an oriented version of £. For an oriented edge e, let Ng(e) be 
the number of times C visits e. One has 

N c {e) = N £ {e) + N £ {-e ), 

and, since the distribution of C is invariant under reversal of all loops, 

N E {e) = N E {-tt). 

Hence (Nc(e)) = 2(Ng(e)). 

Fix an oriented edge e. For \t\ < 1, let 


IAi]ei,e2 


•^ei (^L{ei=el T l{ei^e}) if ^ei ^e2 Rlld 1 g] 7^ ^&21 
0 otherwise, 









20 


FEDERICO CAMIA AND MARCIN LIS 


and let Z t denote the partition function of the corresponding loop soup. 


Using expression (3.2) for the partition function and Lemma 4.1 
that 


we see 


= |det~5(Id- A f )| t=1 = 1 |det(Id-A t )| t=1 

det“5 (Id — A) 2 det(Id-A) 

= \ T* [(M - A)- 1 ^A t ] | t=1 = \ Tr [(Id - A)- 1 ^] 

00 1 

= 2 Tr = 2 Tl ' [ J ^ (Id " A)_1 “ Id )] 
n=l 

= \(Ge : e-l), 

where the fourth identity follows from Jacobi’s formula for the derivative of 
a determinant, and where [Ig]g lt g 2 = 1{ g 1= g 2 =g\■ □ 


As in the case of the partition function, using the above result which 
relates the one-point function to the underlying matrix, exact computations 
can be made for the homogenous model. 

Lemma 5.6. The eigenvalues of A^(p) are ±x with multiplicity d— 1, and 


2d - 1 


x 


ip ± y/ a 2 - 2d + 1 


with multiplicity 1, where a p = Yli=i cos (ifP*)- 

Proof. The result follows from Lemma |5.1| □ 


Note that it follows that the spectral radius of A d is equal to (2d — T)x 
and is achieved by one of the multiplicity-one eigenvalues of A^(p) for p = 
( 0 ,..., 0 ). 

Corollary 5.7. For any edge e of Q^, 


(Nc(e)) 


i \ - i - x Ei_i cQs (ifp.) 

dn d ^ 1 + ( 2 d- l)i 2 - 2xJ2U cos (^ft) 


+ A(x), 


and hence, in the thermodynamic limit, 


Jim (N c (e)) 


1 - x COS OLi 


d(2n) d J[ 0 ) 2 Tr] d 1 + (2d — T)x 2 — 2x J2i= i cos a i 


da + A(x), 


where A(x) = Ap(l — x 2 ) 1 — 1 is smooth on ( 0 , 1 ). 
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Proof. Let cr^(p) be the spectrum of A^(p). Using Lemmas 
have 


5.5 


and 


5.6 


we 


(N c (e)) + 1 = (±Y1 N ^)) + 1 = ^ Tr(Id - A ^) -1 


\\E\ z 

1 1 e'S E 


1 


2\E\ 


Lee a -*)- 1 


2 dn d 


p £Gn Ae<r^(p) 


1 — xa T 


d LE L (i _ x 2 r i , JL V _‘_. 

d dn d 1 + (2d — l)x 2 — 2xa p 


□ 


Corollary 5.8. As x /* x c = 1 /( 2 d— 1), \\m.gd^d(Nc(e)) stays bounded 
for d> 3, and diverges logarithmically for d = 2. 


Proof. Corollary 5.7 and a computation analogous to the one in the proof 
of Corollary 5.4 yield 


rVd 


lim (Nc(e)) ~ 2 d I [p(x) + r 2 ] 1 r d 1 dr. 


Qt/K 

The integral is convergent for all d > 3; for d = 2, one has 
' , v / 2 r( } r ^ fp(x) + 2' 


Jo \p( x ) + r2 \ 

which diverges logarithmically as x /* x. 


= 2 log 


p(x) 


□ 


6. The distribution of Nc(e) 

In this section we compute the probability generating function of Ng(e) 
and then use the result to prove a limit theorem for the two-dimensional 
edge-occupation field. 

For \z\ < 1, let 

OO 

Pe(z) = ^P(A(c(e) = n)z n 

n=0 

be the probability generating function of Nc(e). For a directed edge e, let 
Fg be the partition function of directed loops uj rooted at e which do not 
visit — e and such that = e only for i = 1 and i = |tJ| + 1 (that is, the 
sum over all such loops of the weights of the loops). Let FL be the partition 
function of walks uj rooted at e such that w* = e only for i = 1 and u>i = —e 
only for i = |cJ| + 1 . 


Lemma 6.1. 

rv- zF_g- z 2 flf'_ s \-l* 

V 1 F—g FLF' g ) 


Proof. Let Z z be the partition function of the soup of oriented loops with 
intensity measure ^j2 z , where Jl z (l) = p(l)z Nl ^ e \ One has 
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Using the identity 

I: e£t t e&l t.-eel, e$JL 

and a well-known fact (see, for instance, Lemma 9.3.2 in [l5 or Lemma 4 
in 1221 for a proof in the non-backtracking case), one gets 

(6.1) exp Q MV) \ = (! - G~r l/2 (l - zF _ s )~ 1/2 , 

£: eei 


where Ggg is the partition function of oriented loops w rooted at e such 
that uji = e only for i = 1 and i = |tJ| + 1 , weighted with the weight 
x z (uj) = x(u)z Na<<e \ Splitting these loops according to their first and last 
visit to —e, one has 


Gl s = zF'g{l-zF_g)- 1 zF'_ 


which together with ( 6 . 1 ) finishes the proof. 


□ 


Contrary to dimension three and higher, in two dimensions the edge- 
occupation field is not defined at the critical point (see, for example, Corol¬ 
lary 5.8). Nevertheless, in for any n, including n = oo, one can use 


Theorem 6.1 and the next corollary to prove a limit theorem, as i 1/3, 
for the field normalized by its expectation. 


Corollary 6.2. 


(Nc(e)) = - 


1 F_g+2F'F' 


2 1 - F_g - FLF' - 

e e —e 


Proof. We use the fact that {Nc{e)) = fjfPe(z)|„ =1 


and Lemma 



□ 


Theorem 6.3. Fix n (possibly n = oo) and consider the loop soup C in Q'(. 
Then, for any edge e, as x /* 1/3, Nc{e) / (N (c(e)) converges in distribution 
to the square of the standard normal distribution. 


Proof. We will use Levy’s continuity theorem. Let ip(t) be the characteristic 
function of Nc(e) / (Nc(e)), and let e = 1— Fg— FfF' g&nd C = F_g+‘lF(F'_g. 
From Corollaries 5.7 and 6.2, we can deduce that, asx/’'l/3,£—>-0 and 

by Lemma 


F_g, Fb and F'_g remain bounded. Hence 


6.1 


lim ip(t) = lim 
x/*l/3 e—>o 


1 - 


and Corollary 

(e 2it£ / c - 1 )F_g+ ( e 4it£ /° - 1JFLFWX - l / 2 


6.2 


= (1 - 2 it)~ 1/2 , 


which is the characteristic function of the square of the standard normal 
distribution. □ 


7. The two-point function 

In this section we show the existence of a subcritical regime with expo¬ 
nential decay of correlations. We let 

Nc(e) = N c (e) - (N c (e)) 
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be the truncated two-point function. We denote by Ng(e) the number of 
visits of the loop i to e E E, and write e E £ if i visits e at least once. 

Lemma 7.1. For any pair of edges e,g, we have that 

(N c (e)N c (g)) = E N e (e)Ne(gMt)- 

t: e,g&t 

Proof. We will use the obvious identity 

Nc(e)= E (#W(e) + E (#^)^(e), 

L eSt, t e,g£i 

where jfl is the multiplicity of l in C. Using the fact that £ is a Poisson 
point process, this gives 

<■ N c (e)N c (g )) = <. N c (e)N c (g )) - (N c (e)) (N c (g)) 

= (( E E (#wo?))) 

£:e,g££ £:e,g££ 

-( e E (#w<?)) 

e-.e,gG£ t.e,g£l 

E N £ (e)N e (g) ((#£)(#£'))- E N &) N f.'{ 9 ) <#*> {#*) 

£,£': e,gE£,£' £,£': e,gE£,£' 

= E N e (e)N e (g) (((/ft) 2 ) - (#l) 2 ) = E N t (e)N e (g)»(e). □ 

e,g&£ l\ e,gdl 


On a regular lattice where each vertex has 2d nearest neighbors, the 
number of rooted, no n-ba cktracking walks of length k is bounded above 
by (2d — l) k . Lemma 


7.1 


then makes it clear that for x < 1/(2 d — 1) one 
should expect exponential decay of the truncated two-point function, iden¬ 
tifying x < 1/(2d — 1) as the subcritical regime and x = 1/(2 d — 1) as the 
critical point. 

We now express the truncated two-point function in terms of the Green’s 
function, and use the expression to give a proof of exponential decay in the 
subcritical regime. 


Lemma 7.2. For any pair of edges e,g, we have that 

2(Nc(e)Nc(g)) = Gg } gGg £ + G-g^Gg-g + G_g_gG_g_g+ Gg^gG^gg. 

Proof. As for the one-point function, we consider the soup of oriented loops 
C. We have 


(N c (e)N c (g)) = <(. N £ (e ) + N £ (-e)) (iN £ (g) + N £ (-g))). 

Let us compute (N £ (e)N £ (g)) ■ For \s\,\t\ < 1, let Z s g be the partition 
function of the soup of oriented loops with intensity measure \g s ,ti where 


m?) = ms N ^t N ^ 
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with Nj{e) being the number of visits of i to e. 
partition function (see (3.2)), one has 


Using the form of the 


(N?(e)N.(g)) = 


d_d_ 7 

ds dt^ 1 


S,t 


s,t =1 


Zi,i 


§s §t exp (log z s,t)\ s=t=l 

Zi,i 


d d | d d , 

= ^» logZ -* !«.«=■ + ds' ogZ -<Bt' ogZ - t U«-1 

where we used Lemma |5.5[ and where G s is the Green’s function for the 
non-backtracking walk with weight x s (6j) = x(u))s Na ^ . Let Ggg be the 
partition function of oriented loops Q rooted at g such that tDj = g only for 
i = 1 and i = |tD| + 1 (that is, the sum over all such loops of the weights 
of the loops). Splitting loops that visit g multiple times into loops that 
visit g only once, one can see that = (1 — sGg^) -1 — 1, and hence 
Ggg = Cs( 1 — sGgfi)- 1 , where C does not depend on s (the additional 
factor s comes from the first visit of the walk to y). Note that only loops 
visiting g survive the differentiation Jj, and that 


k s( l - = (1 - g 9 -, 9 -)- 2 = g| 9 . 

Decomposing the rooted loops starting at e according to their hrst and last 
visit to g and using the identity above we obtain that ^G|-| s=1 = Gg^G^g, 
which finishes the proof. □ 


Corollary 7.3 (Exponential decay). For the homogeneous model on the 
torus with x e = x < 1/(2 d — 1), one has 


(N c (e)N c (g )) < 2 


(x(2 d - 1 )) 2d ^a) 
(1 — x(2d — l)) 2 ’ 


where d(e,g) is the graph distance between e and g. 


Proof. By Lemma 7.2 

induced L 1 norm of the transition matrix is 


it is enough to prove that Ggg < • The 

i = x(2d — 1). Hence, 


Ge,j = [(Id - A d n )-% 3 = [(AjI^tld - Af.)- 1 ]^ 

< II(A„) d ( e,9 ((Id - Aj)- 1 !!, < ||Aj||f e ’ 9) (l - IIAJIIO 


-1 


□ 


8 . A REFLECTION POSITIVE SPIN MODEL 

In this section we consider the loop soup on the square lattice, defined 
as the graph with vertices Z 2 and edges between nearest-neighbor vertices. 
The dual graph (Z 2 )* is again a square lattice. To each vertex v* of the 
dual graph (Z 2 )*, we associate a (±1)-valued (spin) variable a v *. 

We define a spin model on the vertices of (Z 2 )* by taking a loop soup C 
in Z 2 and assigning, to each dual vertex v*, spin 

v lec ' 


( 8 . 1 ) 


a v * = exp 
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where 0((v*) is the winding angle (a multiple of ±27 t) of loop t around v* 
(and i here denotes the imaginary unit). Here we assume that there are 
only finitely many loops surrounding any dual vertex. This is true for any 
x < 1/3. 

We will now show that this spin model is reflection positive. (See |5 
for more information on the concept and use of reflection positivity in the 
context of lattice spin models. The question of reflexion positivity in the loop 
soup context is addressed in Chapter 9 of 19 .) Z 2 has a natural reflection 


symmetry along any line l going through a set of dual vertices. Such a line 
splits Z 2 in two halves, Z+ and Z 2 . We also split accordingly the dual 
graph (Z 2 )* in two halves, (Z+)* and (Z 2 )*, such that (Z^_)*n(Z 2 )* = V*, 
where V)* is the set of vertices of (Z 2 )* that lie on l. 

Let F + (respectively, F ~) denote the set of all functions of the spin vari- 


(respectively, (cr v 


_)*; 


Let i? denote the reflection 


ables (av)u*e(z 2 _)* 

operator, d : F± —> F ^, whose action on spins is given by $(a v *) = <r$r v *\, 
and let E denote expectation with respect to the loop soup. 

Lemma 8.1 (Reflection positivity). For all functions f,g G F + , E (fdg) = 
E (gtif) and E (fdf) > 0 . 

Proof. Take v* 6 (Z±)* and draw a path in (Z±)* from v* to infinity. Let 
n v * be the number of loops and arcs of the soup in Z^_ which cross the path 
an odd number of times. One can see that cr„* = (—l) n ”*, and hence the 
spin model in (Z^_)* is a function of the soup restricted to Z^_. 


Let Ei denote the edges of Z 2 crossed by the reflection line l. Theorem 3.3 
implies therefore that / and fig are independent conditional on Ns\e x - Us¬ 
ing reflection symmetry, this gives E(f&g\N s \ E i) = E(f\Ns\ El )E{'dg\Ns\ El )- 
The proof of the lemma follows immediately from the law of total expecta¬ 
tion. □ 

Defining the spin model on the square lattice is convenient to express 
reflection positivity, but one can of course define the same model on other 
two-dimensional graphs. Consider, for example, the spin model defined via 
( | 8 . 1 [ ) on a finite subset of the square or hexagonal lattice with zero boundary 
condition £ and equal edge-weights x e = x Ve E E. We define the (unnor- 
malized) magnetization field as a v*^v *, where S v * denotes the Dirac 
delta at v*. We conjecture that, if x = 1/3 in the case of the square lattice 
or x = 1/2 in the case of the hexagonal lattice, the magnetization field, when 
properly rescaled, has a continuum scaling limit which gives rise to one of 
the conformal fields discussed in Jr] and ( 8 j, namely, the winding held with 
A = 1/2 and f3 = tt in the language of |7J. (Note that, if x < 1/3 on the 
square lattice or x < 1/2 on the hexagonal lattice, the truncated two-point 
function can be easily shown to decays exponentially by arguments similar 


to those in the proof of Lemma 7.1 ) 

Some evidence in favor of this conjecture comes from Q. In that paper, 
the authors consider the scaling limit of a quantity that can be interpreted 


as the one-point function of a spin model defined via ( 8 . 1 ), but using walks 


without the non-backtracking condition. They express the limit in terms of 
the Brownian loop soup |1T] in a way that is consistent with our conjecture. 
We expect that the non-backtracking condition does not change the scaling 
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limit of the loop soup, so the same computation should apply to the one- 
point function of our spin model. 

Our conjecture and Lemma [8TT| suggest that the winding field of (7| and [8| 
with A = 1/2 and (3 = it is reflection positive and satisfies the Osterwalder- 
Schrader axioms of Euclidean field theory [|23|. 

We also conjecture that, both in the case of the square lattice with x = 1/3 
and in the case of the hexagonal lattice with x = 1 / 2 , in the continuum 
scaling limit, the collection of boundaries of the spin clusters converges to 
CLE 4 , the Conformal Loop Ensemble with parameter k = 4. 
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